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Dendronized polymers consist of an elastic backbone with a set of iterated branch struc- 
tures (dendrimers) attached at every base point of the backbone. The conformations of such 
molecules depend on the elastic deformation of the backbone and the branches, as well as on 
nonlocal {e.g., electrostatic, or Lennard- Jones) interactions between the elementary molecular 
units comprising the dendrimers and/or backbone. We develop a geometrically exact theory 
for the dynamics of such polymers, taking into account both local (elastic) and nonlocal in- 
teractions. The theory is based on applying symmetry reduction of Hamilton's principle for 
a Lagrangian defined on the tangent bundle of iterated semidirect products of the rotation 
rt groups that represent the relative orientations of the dendritic branches of the polymer. The 

'^ resulting symmetry-reduced equations of motion are written in conservative form. 
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1 Physical and mathematical motivations 

Dendronized (dendritic) polymers represent an exciting recent advance in the field of chem- 
ical synthesis. These compound molecular structures are formed by assembling multiple 
dendrimers (a low molecular weight unit to which a number of dendrons, or branches, is 
attached) that are each connected by its base to a long polymeric backbone. Known as 
rod-shaped dendrimers, the first patent for synthesizing them was filed in 1987 [ ]. The de- 
scription of dendronized polymers appeared for the first time in the scientific literature in 
1990 [ ]. Since then, the field of dendronized polymers has grown explosively, as a great 
number of papers have been published on the subject. For discussions of recent progress, 
literature reviews, and methodology we refer to [■'», 4, 5]. 

Although the literature describing experimental approaches (such as fabrication and mea- 
surement) of dendronized polymers is extensive, much less work has been devoted to modeling 
their dynamical properties. In our opinion, this dearth of modeling is probably due to the 
mathematical difficulties posed by the analysis of simultaneous forces caused by the elastic 
deformations of the polymer backbone and its attached dendrimers, as well as long-range 
interactions between the dendrimers through screened electrostatic and other forces. From 
the mathematical point of view, obstacles arise in the extension of the classical and well- 
developed Kirchhoff theory of elastic rods [ ], mainly because the description of dendronized 
polymers is geometrically considerably more complex and the long-range interactions among 
different parts of the polymers. We shall note, however, that there has been some recent 
progress in modeling dendronized polymer conformations. For example, the course-grained 
approach [7] and the atomistic/ Janus chain model [ ] have been used to study the formation 
of helical structures. In addition, Monte-Carlo simulations were used to predict stiffness of 
polymer chains in [v]. Nevertheless, we are aware of no previous work that is capable of de- 
scribing the spatio-temporal evolution of the dendronized polymer. This evolutionary aspect 
is important, for example, in the computation of the propagation of sound waves along the 
polymer, or in choosing a dynamical route toward a final configuration. 

This paper aims to derive the general and geometrically exact theory for dynamical evo- 
lution of dendronized molecules from fundamental physical principles. Our theory can deal 
with arbitrary elastic and nonlocal interactions arising from charge distributions on branched 
structures and is thus applicable to a wide range of dendronized polymers. The system we 
have in mind is illustrated on Figure 1.1: rigid bouquets of charges are attached to each point 
of the deformable elastic centerline. In addition, another set of rigid bouquets is attached to 
each end of the first bouquet. We derive equations of motion in a general geometric setting 
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for an arbitrary number of sequential bouquets. However, when transferring these equations 
from their general geometric form to the "nuts and bolts" vector equations for angular and 
linear momenta of each bouquet, these equations become algebraically quite complicated if 
the number of charge bouquets in the sequence is greater than two. 

The description of such a compound system will be made possible by using iterated semidi- 
rect products of Lie groups. The equations of motion are then derived by using a Lagrangian 
reduction by symmetry inspired by the standard Euler-Poincare reduction on Lie groups [10]. 




Figure 1.1: Motivating physical system for the paper: a sketch of dendronized poly- 
mer. Rigid conformations of charges (red) are distributed along the centerline (solid 
black curve). At the end of each branch in a given rigid bouquet, another rigid bouquet 
of charges is attached (denoted here by two blue spheres on a rigid green rod). The 
orientations of the bouquets attached at the ends of different branches of the same 
dendrite (red bouquet) may, in general, be different. More bouquets can then be at- 
tached to the blue spheres, and so forth. The system experiences elastic forces induced 
by deformations of the centerline and the relative deformation of bouquets, as well as 
electrostatic and other nonlocal forces (e.g., Lennard- Jones) due to the interactions 
between each pair of charges. 



Euler-Poincare reduction. Euler-Poincare reduction is a useful tool for deriving re- 
duced Euler-Lagrange equations when a variational principle admits a Lie group of symme- 
tries. In the simplest case, one starts with a Lagrangian L : TG — )• M defined on the tangent 
bundle TG of a Lie group G and assumes that it is invariant under translation by G acting from 
the left, for example. This Lie symmetry allows one to define a reduced Lagrangian £ : g — )• M 
on the Lie algebra g = TG/G of G. Hamilton's variational principle 6 f^ ^ L{g, g)dt = for the 
left G-invariant Lagrangian L{g,g) = L{hg, hg) for any h £ G implies a variational principle 
S f^ ^ £{^)dt = with i{^) := L(e, g~^g) for the reduced Lagrangian, provided the variations of 
^ := g'^^g € g are of the form 5^ = 7]+ [^, -q] for some curve ?? E g satisfying r/(to) = r/(ti) = 0. 
From this equivalence, one shows that the Euler-Lagrange equations are equivalent to the 
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Euler-Poincare equations 

d 6£ _ 51 

where S, = g^^g- This formulation is useful for the evolutionary description of rigid bodies 
and incompressible flows, for example. In [ui] and [i j], this approach has been generalized to 
include the case of systems defined on semidirect products of Lie groups with cocycles in order 
to gain geometric insight on advected quantities. Systems treated by this theory include heavy 
tops, compressible fluids, complex fluids, and liquid crystals, for example. As we will show, 
the description of the dynamics of multi-bouquets summons a more sophisticated version of 
Euler-Poincare reduction involving iterated semidirect products. 

2 Review of the geometry of molecular strands 

In this section, we briefly recall the theory of affine Euler-Poincare reduction and its ap- 
plication to the dynamics of molecular strands that do not have a dendritic structure. Our 
exposition follows the original exact geometric rod theory [12] that has been adapted to molec- 
ular strands with nonlocal interactions in [13, 14]. 

2.1 AfRne Euler-Poincare reduction 

Let G be a Lie group. For the molecular strand G = SE{3), the special Euclidean group of 
orientation preserving rotations and translations in three-dimensional space. Let F be a left 
representation space for G and let c : G — )■ y* be a group one-cocycle. That is, c satisfies 

c{gh) = gc{h) + c{g) , 

for all g,h £ G. Suppose a £ V* and denote by 

0g{a) ■■= ga + c{g) 

the induced affine action of G on V* . Recall that ga G V* is defined by {ga,v) = (^a,g~^v'), 
for all V £ V. 

Consider a Lagrangian La„ : TG — )■ M depending on the parameter oq G V* such that the 
function 

L:TGxV*^R, L{g,g,ao) := La,ig,g) 

is G-invariant. Then i : q x V* ^ M given by 

£{^,eg{ao)):=L{g~^^,ao) 
is well defined. Hamilton's variational principle 

^ / Lao{g,g)dt = 
■Jto 

reduces then to the constrained variational principle 

5 f'l{ta)dt = 0, 

Jto 

for special variations S^ £ Q and 6a £ V* that are obtained as follows. If we denote ^ = 
g~^g £ 01 then a direct computation shows that 

fiTi fin 

s^ = S{g~^g) = -g^^^g g'^g + g'^^g =■ ^ + [C , v] = "^ + ^"^^ ^ > 
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where ry = g^^6g G g and [ • , • ] is the Lie bracket in g, the Lie algebra of G. The variation in 
a £ V* may be obtained from the path a(e) := 0g-i(e)ao G V* , where g{e) G G is a path with 
5(0) = e and g'{0) = r] £ q. Then 

6a = 6(^g^ a + c{g^ )) = —rja — dc(r/) , 

where dc is the derivative of the one-cocycle c at the identity, and r] = g^^dg G g is an 
arbitrary curve on [toi^i] such that ??(to) = ^(^1) = 0. This symmetry-reduced variational 
principle leads to the affine Euler-Poincare equations 



WtTi-^Ti^Ta'"'-^' \Ta, ^^_^^ 

—a + ^a + dciO = 0, 

see [11]. Here o : V x V* ^ q* is defined by {uob,C) ■= - {Cb,u) := {b,Cu) for all (eg, 

u E V, b G V*, where C^ S g denotes the infinitesimal g-representation on V induced by the 

given G-representation. 

The last two terms in the first equation of the system (2.1) come from variations in a as 

follows: 

/ 61 ^ \ I 61 , , A 1 61 . T ( 61 



For the molecular strand, a generalization is needed, namely, we suppose that the La- 
grangian Lqq is known only for a fixed a^ G V* . In particular, we do not suppose that this 
Lagrangian comes from a G-invariant function L : TG x 1/* — )• M. In this case, we assume that 
Lag : TG — )• M is only invariant under the isotropy subgroup G^^^ := {g G G | Og{ao) = a^} of 
ao, that is, 

Lao{hg,hg) = Lao{g,g) 

for all h G G^^. The reduction goes as before, ({S,,Og{aQ)) = Laa{g~^$), except that now 
the reduced Lagrangian d. is only defined on the submanifold g x O^ C g x y*, where 
^ao •" {^g(^o) I g G G} C T^* is the G-orbit of oq. The associated quotient map is, as before, 

(g, g)eTG^ {g-^g, c{g-^)) e g x V* . (2.2) 

The reduced equations of motion are still given by the affine Euler-Poincare equations (2.1), 
where i is arbitrarily extended to a smooth function on q x V* . It is shown in [ ] that these 
equations are independent of the extension of £ to g x y * . Note also that the solution of the 
second equation in (2.1) is a{t) = 0^(()-iOo, where ^(t) = g{t)~^g{t), g{0) = e. 

In many applications oq = 0. It may also happen that the reduced Lagrangian i cannot be 
explicitly expressed as a function of ^ and a only and it has the form £{£,, a, g), where g £ G 
is such that 9g-iao = a. If we assume that i is G^^ -invariant, that is, i{^,a,hg) = £{^,a,g) 
for all h G G'^^, then i{^,a,g) is a well defined function of ^ and a and equations (2.1) are 
still valid, where one computes the functional derivatives as if g were expressed explicitly in 
terms of ^ and a. Even if this is not possible, the derivatives still have explicit expressions in 
many applications, as shown in [l-l] and discussed below. 

2.2 Molecular strand dynamics 

We briefly review the geometric setting for the molecular strand from [13, 14]. Denote / := 
[a, 6] C M and let SE{3) be the special Euclidean group of M^ consisting of orientation 
preserving rotations and translations. Denote by se(3) the Lie algebra of SE(3). The material 
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configuration variables are the spatial position r{s,t) G M^ of the filament and the rotation 
A(s, t) of the rigid charge conformation at parameter value s and time t. The configuration 
space of this system is therefore the group J-{I,SE(3)) of S'£'(3)-valued smooth mappings 
on /. The multiplication in J-'{I, SE(3)) is given by pointwise multiplication in SE{3): 

(Al, ri)(A2, rs) = (A1A2, n + Airs). (2.3) 

The time and space derivatives yield, respectively, the material velocity {A{s,t),r{s,t)) 
and the angular and linear deformation gradients (A'(s, t),r'{s, t)). Given A and r, we define 
the following reduced convective variables [12], 



(2.4) 



The physical interpretation of the variables (2.4) is as follows. The variable p{s, t) represents 
the position of the filament in space as viewed by an observer who rotates with the rigid charge 
conformation at (s,t). The variables ($7(s, t), r(s, t)) describe the deformation gradients as 
viewed by an observer who rotates with the with the rigid charge conformation. The variables 
(a;(s, t),7(s, t)) describe the body angular velocity and the linear velocity as viewed by an 
observer who rotates with the rigid charge conformation. 

Remarkably, the convective variables (2.4) can be obtained by considering the quotient 
map 

TT{I,SE{3)) ^ J'(/,se(3)) x J^(/,so(3)) xT{I,R^f, 

(A,A,r,f) ^ (a;,7,^,r,p) = ((A,r)-i(A, A,r,f), c((A,r)-i)) , (2.5) 

where (A, r)~^(A, A, r, f) denotes the tangent lift of left translation on SE{3) and c is the 
group one-cocycle defined by 

c{A,r):={iA,r)UA,rr\-r). 

This is a cocycle with respect to the representation of the group T{I, SE{3)) on the vector 
space J"(/,se(3)) x ^{1,]^^) given by 

(A,r)(0,r,p) = (Ad(A,,)([7,r),Ap), 

where (A, r) eT{I,SE{3)), {n,T) G -F(/,5e(3)) = -F(/,so(3)) xT{I,R^), p £ T{I,R^), and 
Ad(A,r) is the adjoint action of the element (A,r) G T{I, SE{3)) on se(3). Hence the affine 
representation oi J'{I,SE{3)) on J"(/,sc(3)) x J"(/, M'^) is given by 



7(A,r) 



(f^,r,p) = (A,r)(f7,r,p) + c(A,r). 



It is easy to see that the isotropy group of (0, 0, 0) relative to this affine representation is 
J-"(/, SE{3))q = SO{3). The projection (2.5) is therefore of the general form described in the 
quotient map (2.2), with G = F{I, SE{3)) and V* = J"(/,se(3) x F{I,R?). 
We consider reduced Lagrangians of the form 

£(w, 7, S7, r, p) = ^ioc{uj, 7, ^, r, p) + inp{^, r, p, (A, r)). 
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where the first Lagrangian lioc is explicitly given in terms of the variables {uj,^,^,T, p) 
and the second Lagrangian £„p has still a dependence on (A,r), where (A, r) are such that 
(A~^A', A~^r', A~-'^r) = (f],r,p). The Lagrangian i^p is well defined because it is required 
to be 50(3)-invariant in the following sense. Regard 50(3) as a subgroup of T{I,SE{3)), 
A G S0{3) I— )• (A, 0) E J^{I, SE(3)). Then it is assumed that £np is invariant relative to the 
multiplication (2.3), that is, 

£„p(f],r,p,(/iA,/ir)) =£„p(0,r,p, (A,r)), for all h £ SO{3). (2.6) 

For the molecular strand, this dependence on (A, r) is given in terms of the S'0(3)-invariant 
variables (^, k) defined by 

(as, s'),Kis, s')) := {Ais),ris))-\A{s'),ris')) = (A(s)-1A(/), A(s)-i(r(/) - r(s))) 

and we have the concrete expression 

,{n, r, p, (A, r)) = II C/(e(s, s'), k{s, s'),T{s), T{s'))dsds' 



^np\ 



for a given function U : 5o(3) x (M'^) — )• M. In this expression, the relation (2.6) follows by 
S'0(3)-invariance of the variables ^(s, s') G so(3), k(s, s') G M^. 

Let us show that this general form allows us to treat a nonlocal potential energy of inter- 
action depending on the spatial distances between the individual charged units in different 
locations along the strand. The spatial reference state for the fcth charge in a given rigid 
charge conformation is the sum r{s) + ^/c(s), where r/kis) is a vector of constant length. In 
the current configuration, the position Cfe(s) of the charge k is thus Ck{s) = r{s) + A{s)r]f,{s) 
and the distance from charge k at spatial position Cfc(s) to charge m at position Cm{s') is 

dkm{s,s') := \ck{s) -Cm{s')\. 

The chain of equalities 

dkniis, s') : = \ck{s) - Cm{s')\ = \r{s) + A(s)r/fc(s) - r{s') - A{s')r]^{s')\ 
= \A{sr\r{s) - r{s')) + Vk{s) - A{s)''A{s')Vm{s')\ 
= I - k{s, s') + T7fc(s) - C(s, s')t]^{s')\, 

shows that dkmis, s') can be expressed in terms of the quantities ^(s, s') G so(3), k(s, s') £ M^, 
as required. 

The equations of motion associated with the reduced Lagrangian i are therefore obtained 
by affine Euler-Poincare reduction for a S'0(3)-invariant Lagrangian defined on TT{I, SE{3)) 
and with respect to a reference configuration given by oq = 0. Using (2.1) produces the 
equations, 

, Si ^^ ^ ^ 51 5i ^ 6i 61 

{dt + u;x) ^ + {ds + nx) —- + p X — +rx— + 7X — = 0, 

where u is the vector in M^ defined by the condition ww = cj x w, for any w G M^, and 
similarly for fi. In the equations above, the functional derivatives of ^ relative to (fi, F, p) are 
taken as if one could have expressed explicitly i only in terms of the variables (a;, 7, il, F, p). 
Of course, this dependence cannot be explicitly written down since there is no concrete 
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expression of (A, r) in terms of {ft, T, p) but, as shown in [ ], the derivatives (5£/(5fi, 5i/6T, 
bijbp have an exphcit expression involving only the variables (0^,7, ri,r, p). Moreover, the 
system (2.7) can be written more explicitly as 



(dt + ux) — h [ds + rix) — — - +7 X — hr X ^ — — — ^ + p X 



loc 



6u ^ ' 6fl ' 6f ST "^ 6p 

{dt + ujx) ^— + {ds + nx) — — ^ —- 

07 oi op 

where the term Z{s,s') is the vector given by 

Zis, s') = as, s') (^{s, s')) ^ -^{s, s')eis, s'), 



and the functional derivatives of £np with respect to the variables (fi, F, p) are computed in 
the usual manner. For all the details leading to this system see [14]. 

3 Reduction for iterated semidirect products 

In order to extend the theory of molecular strand dynamics to dendronized polymers, we 
need to utilize additional mathematical concepts, based on the notion of iterated semidirect 
products. Physically, the rotation of the kth bouquet level induces the motion of all subse- 
quent bouquets starting with the level k + 1; this can be formulated mathematically in terms 
of iterated semidirect products. 

3.1 Semidirect product with cocycle 

Let Gi and G2 be Lie groups and assume that Gi acts on G2 by Lie group homomorphisms. 
Define the semidirect product Gi (S) G2 with multiplication 

{91, 92) {91,92) ■■= {9191,92(91 -92)), 91, 91 G Gi,g2,g2 G G2, 

where gi ■ ^2 denotes the action of gi on ^2- The identity element is (e,e) and (51,52)""^ = 
{9i ,9i '92 )• ^s ^^ shall see later, the choice Gi = SE{3), G2 = 50(3), where the 
semidirect product 5£'(3)(DS'0(3) is defined by the action (Ai,r) • A2 := AiA2Aj"\ (Ai,r) G 
SE{3), A2 G SO{3), will be of great interest in the geometric study of two-level bouquets. 
We suppose that Gi (S) G2 acts on a dual vector space V* by an affine action 

^(91,92)0 = {9i,92)a + c{gi,g2), a £ V* 

where a 1— )• {gi,g2)a denotes a representation of Gi(S)G2 on V* and c : Gi(DG2 — )■ F* is a 
group one-cocycle with respect to the above representation, that is, 

c{{gi,g2){gi,g2)) = (51,52X51,52) + £(51,52) 

for aU (51,52), (51,52) G Gi®G2. 
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3.2 Reduced equations of motion 

Let L : r(Gi(s)G2) x F* — ?• M be a Gi-invariant Lagrangian under the affine action of Gi on 
T{Gi(§)G2) X V* given by 

{91, 91, 92, 92, ao) ^ {hgi,hgi,h- g2,h- g2,{h,e)ao + c{h,e)), h G Gi. 

The quotient space (r(Gi(DG2) x V*) jGx is in a natural way isomorphic to the space gi x 
TG2 X y*, the associated projection being given by 

(51, 51, 52, 92, ao) e T(Gi(DG2) x y* 

^ (fi'r^5'i,9r^-5'2,5'r^ •52,6'(gi,e)-iao) Gfli xrG2 x V* . (3.1) 

This reduction map is a special case of the metamorphosis reduction described in [ ] , with the 
exception that here one has an additional parameter oq that produces an advection equation 
in the convective description. 

In our particular case of two Lie groups, there is an additional isomorphism that identifies 
the quotient space gi x TG2 x V* with a simpler space, namely 

31 X TG2 X y* ^ 01 X G2 X 02 X y*, {bJx,V,V,a) ^ (wi,P,P~^P,6'(e,p)-ia) . 
By composing this isomorphism with the quotient map (3.1), we get the projection 
(91, 91, 92, 92, ao) ^ {91^91, 9i^ -92,91^ ■ 92,0(g^,e)'^aQ) 

^ (91^91, 9i^ ■92,{9i^ ■92y^{9i^ •52),6l^^^^_l^^^j-l6'(g^,e)-lao 

= {9i^9i,9i^ -92,91^ ■ (92^92), O^g^^g^)-iao) (3.2) 

since (51, e) (e,5J~ • (72) = {91,92)- This quotient map will be shown to produce appropriate 
convective reduced variables in our applications. We denote by 

^ : 01 X G2 X 02 X F* ^ M 

the associated reduced Lagrangian defined by 

L{9i,gi,92,92,ao) = i{9i^9i,9i^ -92,91^ • i92^92),0(g^^g2)-iao) =: i{uji,p,uj2,a), 

where we have introduced the reduced variables defined by 

^1 = 9i^9i 

P = 9i^ ■ 92 
^2 = 9i^ ■ 92^92 

« = ^(9i,92)-i"0 = ^{gi,grp)-iao- 

In order to compute the reduced Euler-Lagrange equations, we calculate the variations of 
these variables. A direct computation yields 

6uji = f]i + [uji,7]i], r]i = g^^5gi 
6p = pr]2 - rji ■ p, r]2= 9i^ ■ g2^5g2 

6UJ2 = ??2 + [W2, r/2] + UJl ■ rj2 - 'qi ■ 002 

5a = -r]ia - r]2a - (9ic(r?i) - d2c{r]2). 
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where r/j are curves in Qi vanishing at the endpoints and diC : Qi — t- V* , i = 1,2, are the 
partial derivatives of c, that is, 



diciCi) :-- 



d_ 
di 



c(exp(t^i),e), and 92c(6) : = 



t=o 



di 



c(e,exp(t^2))- 



t=o 



The dots ( • ) appearing in the above expression denote the operations naturally induced by 
the action of Gi on G2: 



gi -^2 



d 
dt 



gi ■ exp{tu}2), 



UJl ■ 0J2 



t=0 



d_ 
'dt 



exp(tC(Ji) • u}2, 



(3.3) 



t=o 



where (71 G Gi, a;i € gi, u;2 G 02- The formula for 8002 is proved in the following way 
5u}2 = S{g^^ ■ 92^92) = -gi^Sgig^^ ■ g2^g2 - gi^ ■ 92^^g2g2^g2 + dtig^^ ■ g2^^g2) 

+ gi^gigi^ ■ g2^^92 + gi^ ■ g2^92g2^^92 

= -r]i ■ UJ2 + r]2 + ^1 ■ ■112 + [^2, V2]- 

To obtain the formula for 6a we recall that if a Lie group G is represented on V, then for 
a = g^^aQ + c{g~^) G V*, we have 5a = —rja — dc{rj), where rj = g~^6g. In our case 
G = Gi (D G2 so that taking g = {gi,gi • p) we get 

g'^^g = (51,51 •p)"^^(5i,5'i -p) = {gi^,p'^){Sgi,^gi -p + gi-h) 
= {gi^^gi,p'^{gi^ ■ {^gi -p + gi- {pm -vi -p)))) =■ iniim)- 



Since 



(m,^2)a 



dt 



(exp(tr/i),exp(t?72))a 



t=o 



dt 



(e,exp(t772))(exp(tr/i),e)a = r]ia + r]2a, 



t=o 



the formula for 6a follows. 

In order to write the reduced equations of motion, we need to define the following diamond 
operations. The diamonds associated to the representation of Gi (S) G2 on V* are 



Oj : F X F* -> g*, {v Oj a, iji) = - {rua, v) = {a, rnv) . 
The diamond associated to the representation of Gi on g2 is 

012 : g2 X 02 ^ q[, {^2 012 ^2, ^1) = (/W2, Vi • 6) • 
Using the above constrained variations in 

t2 



(3.4) 



6 / £{L0i,p,iAj2,(i)dt = 0, 
Jti 

a long but straightforward computation finally yields the reduced Euler-Lagrange equations, 



( d 6i _ 6e f6i 



6i 6i ^ r (5^ 

^2 012 T h ^ oi a - dic — 

0UJ2 6a \6a 



d 61 



6i 



6i 



ad* -^^ -'—-uo ■— —oa-dJi — 
dt 6uj2 ^^ 60J2 6p 6uj2 6a \ 6a 



(3.5) 



to which we add the two equations 

p = puj2 — coi • p and a + uiia + uj2a + dic{uji) + 820(002) = 0. (3-6) 
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Here J12 : r*G2 — )■ 0^ is the momentum map of the cotangent Hfted action of Gi on G2 given 
by 

(Jl2(ag2),6) = (092 '6 -92) , 

where ag^ e T*^G2 and ^1 G 0i. 

By using the expression of the infinitesimal coadjoint action of the semidirect product 
Gi (D G2 , namely 

ad(^-^,^2)(^i,/i2) = (ad(^-^ /^i - a;2 012 /W2,ada;2 ;U2 - ^i • /i2) , 

for wi G 01, L(J2 G 02, /Ui € 0*, |i2 G 02' equations (3.5) can be written in a compact form as 

d_ f 5i_ _M_ 

Conservative form. Given two curves /i(t) G 0* and g(t) G G, we have the formula 

^ Ad;(,)_i /x(t) = Ad;(,).i (^^(t) - ad|(,) /x(t)) (3.8) 

where ^(t) = g{t)^^g{t) G 0. Suppose that the Lie group G acts on a vector space V by the 
representation v 1— )■ (^rt;. Consider the afiine representation a h^ ga + c{g) of G on y*, where 
o I—)- 5a is the contragredient representation and c a group one-cocycle. Then we have the 
relation 

dc^ {gv) = Ad;_i (dc^ («) - ^ o c(<7-^)) • (3.9) 

By using these formulas and supposing that ao = 0, equations (3.7) can be rewritten in the 
form of a conservation law as 

This is easily seen since a evolves as a = O^^g^g^yia^ = (51,52) ""^oo + c((gi,52)~^)- 

3.3 Hamiltonian formulation and Poisson brackets 

In this paragraph, we briefly describe the Hamiltonian side of the reduction process. 

Let H : T*{Gi(§)G2) x F* -^ M be a Gi-invariant Hamiltonian. One can think of H 
as being obtained from the Lagrangian L by a Legendre transformation, the variable in 
V* being viewed as a parameter. On r*(Gi(DG2) x V* we consider the Poisson structure 
given by the canonical symplectic form on T*(Gi(S)G2) and the trivial Poisson structure 
on V* . In this way, Hamilton's equations for H are equivalent to the canonical Hamilton 
equations for Hag together with the equation clq = 0. As on the Lagrangian side, we define 

Haoioigi,0!g2) := -H'(agi, "ga' «o)- 

As in §3.2, we consider the quotient map 

r*(Gi(DG2) X y* ^ 0* X G2 X 0* X V*, 

given by 

("9i'«92'«o) H> (ni,p,H2,a) := {g^^ag^,g^^ -92,91^ • {92^0ig2)^^{gi,g2)-^i'^o)) (^-H) 
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which defines the reduced Hamiltonian h = /i(/ii,p, /U2,a) on 0^ x G2 x 02 x V^*- H H is 
associated to a Lagrangian L by Legendre transformation, the reduced Hamiltonian h can 
also be obtained from i by the Legendre transformation, namely, 



h{f^i,P,fJ-2,a) = (/Ui,a;i) + {^12, '■^2) - i{uji,p,u}2,a), 



5ui 



Aii- 



(3.12) 



To obtain the reduced Hamilton equations, one first needs to compute the Poisson bracket on 



0^ X G2 X 02 X y* obtained by Poisson reduction (see e.g. 
yields 



16, §10.7]). A direct computation 



{/,M 



Ml, 



5f 6h 



/^2 



6f 6h 6f 6h 
' dni 5a 5fi2 Sa 



6f Sh 

6fi2 ' (5/X2 

6h 6f 6h 5f 
5fii 6a 6fi2 5a 



+ {dic 



5fii 



+ 



SI 
Sin 



,J 



12 



+ 820 

6h 
Sp 



si_ 

<5/i2 



6h 
' 5a 



+ 



5_l 6h_ 
5p ' (5/U2 



die 
6h 



/^2 



6h 
6fii 



SI 
6fii 



6h 

SH2 



6h 

Sm 



SI 

SlJi2 



J 



12 



+ d2C 

5h 
5p 



6h 

SH2 



s_l 

5a 

6h 6j^ 

Sp ' 6^2 



The first line represents the Lie-Poisson bracket of the dual of the semidirect product Lie 
algebra 0i (8)02; when the second line is added, we get the Lie-Poisson bracket on the dual of 
the semidirect product (01 (D02) (D V. The third line is due to the presence of the affine term. 
The last line is a new expression arising from the fact that reduction is carried out for the 
subgroup Gi C Gi(DG2 and not the whole semidirect product. The Hamiltonian equations 
associated to this Poisson bracket are 



d 
di 

d^ 
dt 
d^ 
dt 
d 

-rC 

dt 



m = a.d^h_ fJ.1 + J12 

Sfj.1 



t^2 



p = p 



ad*gft ^2 



5h 

Sf^2 

6h 



p 



6h\ 6h 

Sp) Sfl2 

i6h 6h 



5p 5fj,i 



■ f^2 



6h 
Sm 

_5h^ 

(5/ii 6fi2 



P, 



dic 



6h 
SfJ-i 



012 fJ-2 

5h 
5a 



820 



6h 



oi a + dic 



T 



6a 
02 a + 820^ 



6h 



5h 
6a 



6f\ 
6a r 



(3.13) 



6fl2 



If h is obtained from a Lagrangian £ by the Legendre transformation (3.12), these equations 
can be obtained by substituting the relations 



61^ _ ^ __^ 

6fii ' 6p 6p^ 

into the reduced Euler-Lagrange equations (3.5), (3.6). 



6h 
6a 



6a^ 



3.4 Affine reduction at fixed parameter 

In applications, one often has the situation that the Lagrangian La^ : r(Gi(S)G2) — )■ M is 
known only for a fixed value ao £ V*. In this case, the theory developed so far does not apply, 
since I is not the reduction of a Gi-invariant function on T(Gi (DG2) xV*. It turns out that 



the theory in Section 3.2 is applicable, however, if we assume that L^g 
where 

(GiYao ■= i9i e Gi I 9{g^,e)ao = ao}- 



is (Gi 



-invariant. 
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It also happens in applications (this is the case for the bouquets) that the reduced Lagrangian 
I cannot be explicitly expressed as a function of the variables Ci , p, ^2 ; ^ but it has the form 
^(Ci;P)^2;a,5i)) where gi G Gi is such that c^g^ ,p~^) = a. 

If we assume that £ is (Gi)^p-invariant, that is, i{Ci,p,C2,ci',hgi) = i{^i,p,(,2,a,g) for 
all h G (Gi)Jjjj, then £{^i,p,^2,a,g) is a well defined function of {Ci,p, 6,2,0,). Indeed, if 
c{gi ,p^^) = c{h^ ,P^^), then 0, -1 „-i\0 = ^(^-1 p-i)0 which is equivalent to 

{h^\p-')-\gi\p-') = ihg^\e) G (Gi(s)G2)§, 

that is, higi G (G'i)q. In this case, equations (3.5) are still valid, where one computes 
the functional derivatives as if gi were expressed explicitly in terms of S,i,p,£,2,cl. Although 
this may not actually be possible, the derivatives may still have explicit expressions in many 
applications. We refer for details to [ i]. 

3.5 Generalization to N groups 

In this section we introduce iterated semidirect products that will be crucial in the description 
of the dynamics of multi-bouquets. This is motivated by the fact that, as we have mentioned 
earlier, the rotation of a bouquet at the kth level induces, via actions of semidirect products, 
rotation of all subsequent levels of bouquets, and thus the distance between the charges is 
affected correspondingly. The semidirect product rotation is a mathematical property and 
it is applicable to any dendronized polymer with a tree- like structure of the bouquets. This 
representation is independent of the details of the polymer, its chemical composition, etc., as 
it represents the purely topological aspects of the polymer's geometry. 

Note that, in principle, our construction fails should some of the bouquets close at a 
subsequent stage (so the dendrimes become a graph and not a tree). Such graph-like molecules 
have not been synthesized yet. There are also substantial mathematical challenges in dealing 
with graph-like structures as opposed to tree-like structures. In this paper we will study 
exclusively molecules with tree- like structures. 

Semidirect product structure. We shall present the construction of the configuration 
space of A^- level bouquets, also called A^-bouquets. This will be an iterated semidirect product 
of A" Lie groups modeling the tree structure of a multi-bouquet. 

Let Gi, . . . , Gn be Lie groups. Assume that Gn~i acts on Gn by group homomorphisms 
and form the associated semidirect product group Gn-i®Gn- Then suppose that the group 
Gn-2 acts on Gn-i®Gn by group homomorphisms, that is, 

gN-2 ■ {{9n-i,9n-2){9n-i,9n-2)) = {9N-2 ■ {9n-i,9n-2)) {9N-2 ■ {9n-i,9n-2)) 

and form the semidirect product Gn-2®{Gn-i®Gn)- Proceed inductively and form the 
iterated semidirect product group 

Gi®{G2®{G:,® {■■■ ®Gn) ■■■)) B {91,92, ■■■ ,9n). (3.14) 

For example, suppose that G acts by group homomorphisms on two other groups H and 
K. Then 

g-{h,k):={g-h,g-k), geG, he H, keK, 

is an action oi G on H®K hy group homomorphisms, provided the following compatibility 
condition holds: 

g-{h-k) = {g.h).{g.k). 
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A simple example of such a chain of semidirect products is obtained by considering the 
same group G, acting on itself by conjugation: G(§){G®{G®{--- (S)G) • • • )). 
Returning to the general case, the group multiplication reads 

{91, 92, -, gN){hi, h2, .:, Hn) = {9ihi, 92i9i- h2), g3i92- {91- h)), -, 9n{9n-i- ■■■■ 92- (gi- fiN))) 
and the inverse is 

{9i,92,---,9n)~^ = {91^, 9i^ -92^,91^ • fe^ ■ 9^^), ■■■,9i^ ■ {92^ ■ ■■■■ (Stv-i 'Siv^)))- 

Lagrangian reduction. Suppose that the semidirect product (3.14) acts on the dual V* 
of a vector space V by an affine representation 

OH- 6'(gi,...,gjv)« = i9i,-,9N)a + c{gi,...,gN), a E V* . 

This includes all parameters given by the specific physical situation in the A'^-bouquet. 

Let L : T (Gi (D (G2 (D (G3 (S) (• • • (S) Gn) • • • ))) x ^* -^ I^ be a d-invariant Lagrangian 
under the affine action of Gi on T (d (S) (G2 (D (G3 (S) (• • • (S) Gat) ■■■))) xV* given by 

(51, 51, 52, 52, •••, 5Af, gN, ao) ^ {hgi, hgi, h- g2,h- g2, ..., h- gN,h- g^, O^h,e,..,e)ao) , 

h £ Gi. We consider the quotient map 

T (Gi (s) (G2 (D (Gg© (• • • (DGjv) • • • ))) X y* ^ 01 X G2 X 02 X ... xGnxqn 

given by 

{gi, gi, g2, 92, 93, gs, ■■■, gN, gN, ao) H> {uJi,P2,UJ2,P3,UJ3,---,PN,^N,a) 

■= {gi^gi,9i^ ■ g2,gi^ ■ g2^g2,gi^ -92^ -93,91^ -92^ ' 93^93, ■■■,0{gi,g2,...,g,^)-^ao) . 

We thus obtain the convective variables 

^i=9i^ ■ 92^ ■ ■■■■ gi^gi and Pi = gf^ ' ^2"^ ••••• 5il\ • Ri- 
nsing the semidirect product multiplication, we can write 

{lOi,UJ2,---,^n) = {gi,92,---,9N)~^{gi,g2,---,9N)- 

Note that the group Gi©{...®Gn-i) naturally acts on the group G2®{---®Gn)- The 
action is given by the same expression as the group multiplication on the semidirect product: 

{gi,-,9N-i)* (/i2,-",/iAf) := (51 • ^2,52 • (51 • h:i),...,gN-i ■ {gN-2 ■ ••• ■ gi ■ fiN)) ■ (3.15) 
Using this action, the variables pi can be obtained by the formula 

{p2,-,pn) = {gi,-,gN-iy^ • (52,---,5Af)- 

Reduced Euler-Lagrange equations. A straightforward computation shows that the 
constrained variations are 

d *~^ 

S^i = -Q-rii + [wi, r]i] + ^{uJk -rji-T^k- Wj), r/j := ^f ^ • {g^^ ■ (... • g^^Sgi)) (3.16) 

fc=i 
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and 

i-l 

Spi=Pi'ni-'^r]k-pi. (3.17) 

fc=i 

where rji are curves in Qi vanishing at the endpoints, i = 1, . . . ,N. The dots in the above 
expression denote the operations naturally induced from the action of G^ on Gj, k < i, as in 
(3.3). 

For example, the formula for 6iOi is proved in the following way 

6uji = S{g^^ ■ g^^ ■ .... ■ g^^gi) 



(Si ^^9i9i ^ ■ 92^ ■ ■■■■ ■ 9i ^9i) - (5i ^ • 92 ^^9292 ^ ' ■■■ ' 9i ^9i) 
- ••• - (s-r^ • 92^ ■ ■■■ ■ 9^^^9i9i^9i) + {9i^ ■ 92^ ■ ••• • 9^^^9i) 



i—l r, 

o 



^r]k-uji- rjiUJi + jrAgi^ ■ 92^ ■ ■■■■ 9i ^S9i) + (fi^i ^9i9i ^ ■ 92^ ' ■■■■ ' 9i ^^9i) 



fc=i 



+ (5i ^ • 92 ^9292 ^ ■ ■■■■ 9i ^^9i) + ■■■ + {9i^ ■ 92^ ■ ■■■ ' 9i ^9i9i ^^9i 

i—l n i—l 

fe=i k=l 

= -gl'Hi + [^i^ ^i] + X^('^fc -rii-rik- uji) 
fc=i 

Using these variations in 

/•ti 
5 I i{uJi,P2,UJ2,P3,^3, ■■■,PN,(^N,a)dt = 

Jto 



yields the N equations 

Ail- H* — -V — -^—- V /^ — T 

dt6uJt '^^ Suji ^ ^ buJi * bpi J-^^\ ^ ''^ bujk ''' \6pk 

J — 1 K — 2 + 1 

+ ^o^a-^,c^(^), i = l,...,N, (3.18) 

da \oaJ 

where Oi : V x V* ^ q* is given by (3.4), and ioi i < k 

Oik ■ QkX Qk^ 9*i, (& Oik fj'k,'ni) '■= {^J■k,Vi ■ Ck) , ik G flfc, Hk G 0fc, m e g*, 

and Jj/fc : T*Gk — )• 0* is the momentum map of the cotangent lifted action of Gi on G^- 

Among these A^ equations, note that the first and the last of them are simpler, since one 

of the terms involving the summation is absent. Note also that when i = 1, the term p~ j— 

is absent. (The piS start with p2.) 

Using the infinitesimal coadjoint action of the semidirect product allows these equations 

to be rewritten as 

d f 5i di \ ,, f Si 5i \ S£ , T ( ^^ 

ad/ , , , \ T — : •••: T + -^ o a — dc 



dt\bujx'"' bu^t^j ^'^^'■■■''^'^^\5uji'"'' SunJ da \6a 
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where o : V* xV — ^ [fliCD (02(D(- • • (Dbat) •••)]* is associated to the action of the whole 
semidirect product on V* and 

J:T*{G2®{--- ®Gn)) ^ (01 (D (• • • ®QN~i)r 

is the cotangent bundle momentum map associated to the action of G2 (§)(••• ®Gn) on 
G'i(D(- • • ®Gn-i) defined in (3.15). It is given by 

/ N N \ 

J'(ag2,...,agj^) = ( ^Jifc(agJ,..., ^ Jifc(agJ, ..., JAr_iiv(Qg^) j . 

\fc=2 k=i+l / 

Conservative form. If qq = 0, equations (3.8) and (3.9) may be used to rewrite the 
equations of motion (3.18) equivalently in a conservative form as 

4 Application to multi-bouquets 
4.1 Dynamics of two-bouquets 

In order to apply these ideas to the dendronized polymers, we must parametrize the position 
of each charge. For the two-level polymers presented in Figure 1.1, the position of the base line 
is denoted by r(s, t). Then, at every point s, the rotation of the first level of the rigid bouquet 
(colored in red) is measured by Ai e 50(3). In the exact geometric rod approach [12], Ai 
is measured with respect to a fixed coordinate system. In order to complete the description 
of the system, the rotation of the second level of bouquets (marked green) needs to be taken 
into account. In the exact geometric rod theory, the rotation of the second bouquet level 
A2 G 50(3) would be measured with respect to the same fixed coordinate system, rather 
than with respect to the first bouquet. However, we shall depart from this convention and 
measure the orientation of each successive charge bouquet with respect to a frame attached 
to the previous bouquet. The latter convention will lead to a semidirect product action by 
homomorphisms. 

The rotations of the base bouquet Ai depend on the point on the base curve s and, of 
course, time t. The rotation of the second bouquet level A2 also depends on the numbered 
vertex to which it is attached. We call this vertex ki for s and mi for s' . The position of a 
given charge depends on the number of the charge in the second bouquet (called k2 for s and 
m2 for s'), on the orientations Ai and A2^ and on the position of the base r{s). Thus, A2 is 
measured with respect to a frame attached to the first bouquet and we find the formulas 

Ck.k^is) = r{s) + Ai{s)'nk^{s) + A'^'{s)Ai{s)'n,,^{s), 
Cm,m, {s) = r{s') + Al {s')rj^^ {s') + A™i {s')Ai {s')r,^^ (s'), 

for the position of charges on a two-bouquet anchored along the centerline at positions r{s) 
and r(s'). Note that the orientation Aj refers to the frame in which a charge bouquet at the 
ith level is rigidly positioned and not to the orientations of individual components of that 
charge bouquet. 

If more levels of the bouquets were present, we would need to introduce orientations A3, A4, 
. . . and corresponding labels fcs, k^, . . . at s, and tjis, tu^, ... at s' . As we shall explain below. 
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this notation is rather cumbersome, especially in the case of A^-bouquets. We shall concentrate 
mainly on two-level bouquets here, as this case illustrates most of the mathematical concepts 
and is the necessary first step in applications. Formulas for A^- level dendrimers will also be 
derived and discussed later in enough detail to make the general pattern of the equations 
evident. 

Let us explain how the group SE{3)(§)SO{3) appears already in the study of the two- 
level bouquet problem. We first consider the case of a single multi-bouquet whose kinematic 
description coincides with that of a chain of rigid bodies. 




Figure 4.1: A figure illustrating a chain of coupled rigid bodies connected through a 
ball joint (black spheres). 



Motivation for the semidirect product: coupled rigid bodies. Consider the 
kinematics of two rigid bodies coupled with a hinge joint, as in [ ]. We denote by Bi and B2 
the reference configurations of the bodies. Let Pq be a preferred point (for example the center 
of mass) in the reference configuration Bi. Let Pi be the (fixed) hinge point in the reference 
configuration of the first body and P2 be the (fixed) hinge point in the configuration of the 
second body. Let 77 denote the vector from Pq to Pi in Bi. Fix an inertial frame and denote 
by r(t) the current position of the point Pq. The configuration of the first body is determined 
by an element Ai(t) G 50(3). We denote by A2(t) G S0{3), the orientation of the second 
body relative to the first body. The motion of a point in the first body is thus given by 



x{t)=r{t) + Ai{t)X, X£Bi 



and in the second by 



x(t) = r{t) + Ai{t)ri + A2(t)Ai(t)y, Y G B2. 



(4.1) 



(4.2) 



The configuration space of the system is M^ x S0{3) x SO{3) B (r,Ai,A2). In order to 
show how the semidirect product SE(3)(^SO{3) arises, we consider the kinetic energy of the 
system 



a: Ai,Ai,r,f,A2,A5 



1 



\±{t)\^d-'X + 



1 



1 



Bi 



r + AiX 



' d?X + 



2 j3^ 



B-z 



\mv<^'y 



B2 



r + A1T7 + A2Aiy + AaAiF 



d^y, 
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where d X and d Y denote the mass densities of the corresponding bodies. One observes 
that the kinetic energy K is invariant under the cotangent hft of the SE{3) action given by 

(Al, r, A2) ^ {hAi, hr + v,h- A2) , 

where /i • A is the action by homomorphisms defined by conjugation h ■ A := hAh^^. The 
configuration space is therefore naturahy given by SE{3)®SO{3). By SE{3) invariance, one 
can write the kinetic energy in terms of the reduced variables 

on T{SE{3)®SO{3))/SE{3) ~ se(3) x SO{3) x so(3). If the Lagrangian of the system is 
only S0{3) invariant one has the additional reduced variable p = A^^ r; such a situation 
arises in the case of multi-bouquets. In order to obtain the equations of motion, we apply the 
framework described in Section 3 with Gi = SE{3), G2 = SO{3), and the cocycle is absent. 
In the case of 5-E'(3)-invariance, the abstract system (3.5) gives 

\ 6i 61 51 51 _i ,6£ 

dt + u;ix )- h7 X ^ +a;2 X \- ^p -p ^ =0, 

/ ouji o-f 0U2 op op 

ai + cjixj— = 0, (4.3) 

\ 51 61 .61 

dt + U>2'X T h CJi X p '— =Q. 

/ 00^2 0U^2 op 

The variables used here differ from those used in [ L < ] where the configurations Ai,A2 G 50(3) 
of the bodies are given relative to the same inertial frame. That is, we have the relations 

Al = Al A2 = A2A1 

and (4.1), (4.2) read 

xit) = r{t) + Aiit)X, XeBi resp. xit) = r{t) + Ai{t)T] + A2it)Y, Y e B2. (4.4) 

Note that when (Ai,A2) 1-^ {hAi,hA2h~^) then (^1,^2) 1-^ {hAi,hA2), consistently with 
the symmetry considered in ['• ']. The reduced variables ili = A^ Ai, Q.2 = -^2 ^2) and 
A = A'^ A2 used in [ 1 7] are recovered from our variables wi , 0^2 , P via the simple relations 

r^i = wi , 

O2 = ^^^^2 = (A2Ai)~^ (^A2Ai + A2A1) =0Ji+0J2, 

A = A^^A2 = A5"U2Ai = p. 

On the Hamiltonian side, one has a Poisson diffeomorphism between the two representations. 

Remark 4.1 In [7], the coupled system is actually described slightly differently than in 
equation (4.4). Namely, a point in the first body, Xi ^ Bi, evolves as w{t) + Ai{t)Xi, while 
a point in the second body, X2 G B2-, evolves as if (t) + A2{t)X2^ where w(t) is the current 
position of the hinge between the two bodies. In order to simplify the comparison with two- 
bouquets, we have changed the description slightly by using (4.4) instead, since this approach 
generalizes easily to a chain of N rigid bodies, as we discuss below. ♦ 

In the case of a chain of N coupled rigid bodies, the current position of a point in the ith 
rigid body is 

x{t) = r{t) + Ai(i)T7i + A2(t)Ai(t)r72 + ... + Ai{t)...Ai{t)Xi, Xi G Bu 
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where Bi is the reference configuration of the ith body and iji is the vector between the two 
hinges in the ith rigid body. Therefore, the kinetic energy is given by 



N „ 2 



d-^X: 



i^(Ai,Ai,r,f,A2,A2,...,A7v,A7v) = 77 V / r + Airji + - + dt{Ai...Ai) X, 

and one observes that K is S'£'(3)-invariant under the action 

(Al, r, A2, ..., Aat) ^ (/lAi, hr + v,h- A2, ..., h ■ An) , 

in which the group SE(3) acts by homomorphisms on the other SO{3) groups. The configu- 
ration space is therefore naturaUy given by the iterated semidirect product 

SE{3) (D {S0{3) (D {SO{3) (D (• • • (D50(3)) •••))• 

Remark 4.2 Another difference between coupled rigid bodies and our system of dendrite 
polymers is that each rigid body in the chain (say, at the position i) is completely characterized 
by its orientation matrix Aj. For a dendrite, each branch of the bouquet at each level may 
have a different orientation. For example, suppose that the first level rigid bouquet consists of 
ni branches rjj^^, ki = 1, ..., ni whose orientation at s is given by Ai(s). At each extremity of 
these branches, there is a second level rigid bouquet given by ^2 branches ry^^- The position 
of each of the ends of these branches is then 

Cfci =r{s) + Ai{s)r]k^{s). 
Since these second level rigid bouquets can evolve independently we would need rii variables 
therefore 



A2^(s), ki = 1, ...,ni instead of simply A2(s). The position of a charge at the second level is 



CkMs) = r{s) + Ai{s)rjk^{s) + A^^'{s)Ai{s)rjk^{s). 

Similarly the position of a charge at the i-th. level depends on the branches at lower level. 
Thus if we denote by Uj the number of branches at the jth level and by kj = 1, ...,nj the 
indices of these branches, the variables needed at the ith level are A^'"^' '^^'"'' ^ (^g) instead 
of simply Aj(s). The position of a charge at the i-th level is then 

Cfeite...fc„ = ^(s)+Ai(s)r7fc^(s) + A^iAi(s)T7fc^(s)+ 

.... + Af— ••''=^(.)Aj--'^H«)--Ar(«)Ai(«)^fc.(^) 
The configuration space is then given by 

SE{3) (D (50(3)"i (D (50(3)"i"2 d) (. . . d) 50(3)"i"2-niv) . . . ))^ 

where S0{3)'^ = SO{3) x ... x SO{3), and the product is taken n times. Strictly speaking, 
we should include the variables A-""^^'" ^ in the formulas but this would make the resulting 
expressions appear unnecessarily complex, especially in the case of A^-bouquets. Thus, we 
shall only deal with one variable Aj(s) at the ith level and leave the general case to the 
interested reader. We hope no confusion will arise from having made this simplification for 
the sake of streamlining the notation. ♦ 
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Motivation for the afRne Euler-Poincare theory. Recall that for the molecular 
strand (corresponding to a one-bouquet; see §2.2) the main goal was to treat a nonlocal 
potential depending on the separation distance 

dkm{s,s') ■.= \ckis) -Cmis')\, where Ck{s) = r{s) + A{s)r],.{s). 

The crucial property in the implementation of the affine Euler-Poincare theory is the invari- 
ance of the expression for dkmis, s') under the action of 50(3) given by 

{A,r)^{hA,hr), heS0{3). 

For the generalization to two-bouquets, the position of the charge is given by 

Cfcife2(s) = r{s) + Ai{s)r]k^{s) + A2{s) Ai{s)r] ^^{s) , 

and we need to treat a potential depending on the distances 

"'kik2,mim2 1,^, S j . — IC/j-^fcj \^) f^mim2 \^ j r 

In this case, we have 50(3)-invariance of the distances under the transformation 

(Ai(s),r(s),A2(s)) ^ {hAi{s),hr{s),hA2{s)h-^), h G 50(3), 

which shows that the action by homomorphisms again arises naturally, as in the case of 
the coupled rigid bodies. It is therefore appropriate to consider the semidirect product 
5-E'(3)(S)50(3). On the other hand, for the simpler case of the molecular strand, we al- 
ready know that a cocycle is needed to understand geometrically the reduction process from 
the standard Euler-Lagrange equations to the reduced equations of motion. These observa- 
tions strongly suggest that the geometric setting developed in Section 3 is appropriate for the 
description of multi-bouquets, as we now show. 

Let SE{?>) act on 50(3) by (Ai,r)-A2 := AiA2A^\ where Ai, A2 € 50(3). The semidirect 
product 5-E(3) (8)50(3) associated to this action by homomorphisms of 50(3) has therefore 
the multiplication 

(Ai,r,A2)(Ai,r,A2) = (AiAi,r + Air, A2AiA2A^^) 

and inverse 

(Al, r, A2)-^ = (Ar\ -Ar V, Ar^A2 ^Ai) . 

Formulas for the semidirect product. The conjugation in 5i?(3)(D50(3) is given 
by 

AD(A^^^^A2)(^i'^'^2) = (Ai,r,A2)(Ai,r,A2)(Ai,r,A2)"^ 

= (AiAiA^i, r + Air - AiAiA^V, A2AiA2AiA^1A^1AiAiA^1) . 

Consequently, the adjoint action on of 5ii^(3)(S)50(3) on se(3)(S)so(3) is expressed as 

Ad(A^_^^A2)('^1.7,'^2) 

= (AitJiA^\ Ai7 - AiwiAj; V, A2Ai(a;i + oj2)A'{^A'^^ - AitJiA^^) (4.5) 
and the Lie algebra bracket is given in vector form by 

ad(^^,-),,^2)(cji,7,cj2) = [(0^1,7,^2), ('^1,7,^^2)] 

= (cJi X cJi, 0^1 X 7 — 0^1 X 7, a;2 x 0)2 + ^^i x a>2 + '^2 x '^1) • (4-6) 

Therefore, the infinitesimal coadjoint action then reads 

ad(^i,-y,a;2)(/^i''^'/^2) = (/xi X cji -7 X i/ + /i2 X a;2,i^ X cJi,/X2 x (o^i -h <x)2)) • (4.7) 
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The reduced variables. For the problem of two- level bouquets the configuration space 
J-"(/, SE{3) (§) SO{3)) consists of smooth functions of an interval /with values in SE{3)(§)SO{3). 
This is a group relative to pointwise multiplication and hence we may apply the formulas de- 
rived above. To simplify notation, we use the same symbols (Ai, r, A2) to denote variables in 
J^{I,SE{3)(§)SO{3)). These are functions depending on the curve parameter s £ I. 

The theory developed in Section 3.2 now applies. Take Gi = T{I, SE{3)) 3 (Ai,r), 
G2 = T{I, SO{3)) 3 A2, and the action of Gi on G2 by group homomorphisms given by 
(Ai,r) • A2 := AiAaAj^^ The representation space is F = F* = T{I,R^)^ 3 {fli,T,n2,p) 
where we identify V and V* by the L^-inner product on /. The J"(/, SE{3))(§)J^{I, 50(3))- 
representation on J^{I, M'^)^ is given in vector form by 

(Ai,r,A2)(ni,r,n2,p) : = (Ad(A,,,,A2)("i,r, ^2), Aip) 

= (Ai^i, AiF - Aifii X r, A2Ai(ni + ^2) - Ai^i, Aip) 

and the group one-cocycle c on J-"(/, SE{3))(^J^{I, S0{3)) is taken as 

c(Ai,r,A2):= ((Ai, r, A2)a,(Ai,r, A2)-\ -r) . (4.8) 

The choice of this cocycle is motivated by a convenient expression of the reduced variables, 
as we shall see below. 

Indeed, using the formula of the quotient map T {F{I, SE{3))®T{I, SO{3))) xT{I, M^)^ -> 
J"(/,se(3)) X T{I, 50(3)) x J"(I,so(3)) x T{I,R^)'^ given in (3.2) with ao = 0, we obtain 

(Ai,Ai,r,f,A2,A2) ^ ((Ai,r,A2)"^(Ai,Ai,r,f,A2,A2),A5'U2Ai,c((Ai,r,A2)"^) 
= (wi,7,a;2,P,f^i,r,J72,p), 
where 

(wi,7,a;2) : = (Ai, r, A2)"^(Ai, Ai, r, f , A2, A2) 

= (A-i, -A- V, Ar^A2-Ui)(Ai, Al, r, f, A2, A2) 

= (A^Ui,A^if,A^iA^iA2Ai' 

and 

(f^i,r,f}2,p): = c((Ai,r,A2)-^) = ((Ai,r,A2)-^a,(Ai,r,A2),A^V) 

= (Ar^A;,ArV', Ar^A2-iA'2Ai, A^V) . 

We continue to use the convention that for a vector v € M'^ we denote by f := v G so (3) 
the matrix defined by vw = v x w. Summarizing, we have defined the reduced convective 

variables 

wi = A^Ui Qi = A^U; 

7 = Ar^f r = ArV' 

W2 = A^U^^A2Ai ^2 = A^U^U'a^i ^ ' ' 

p = A^^A2Ai p = Aj^V. 

These variables will be shown to be well adapted to the description of the dynamics of two- 
bouquets. Here, uji are angular rotations of each bouquet, with oji being measured with 
respect to the base coordinate system and uj2 measured with respect to Ai. Similarly, Qi 
and r are the familiar variables describing the infinitesimal twist and stretching of the base, 
whereas ^2 describes the infinitesimal twist of the second bouquet in the coordinate frame 
connected with the first bouquet. Notice that there is no r2 as there is no stretching associated 
with the second bouquet. Finally, the new variable p is the orientation of the second bouquet 
seen from the coordinate frame of the first bouquet. Thus, p is a new variable associated 
with the multi-bouquet structure of the polymer. 
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The Lagrangian for two-bouquets. In this case, the material Lagrangian depends 
only on the variables (Ai, Ai, r, r, A2, A2). It is important to interpret it as being the restric- 
tion of a parameter-dependent Lagrangian defined on TT{I, SE{3)(§)SO{3)) x J^(I, M^)^ for 
the zero value of the parameter in J^(/, M'^)^, that is, 

L(ao=o) =i(Ai,Ai,r,f,A2,A2). 

We are thus in the situation described by the theory developed in §3.4, provided L is {Gi)q- 
invariant. We have (6*1)0 ~ SO{3) C T{I,SE{3)) as a direct verification using (4.8) shows. 

Nonlocal variables. The most general form of the symmetry-reduced Lagrangian for 
level-two bouquets is 

i = £;oc(^i,7,^2,p,fii,r,fi2,p) + 4p((^i,7,<^2,p, J^i,r,ri2,p, (Al,r)), 

where the first Lagrangian ii^c is explicitly given in terms of the variables cji, 7, 0^2, p, ^1, T, 
fl2, p and the second Lagrangian i^p has still a dependence on (Ai, r), where (Ai, r) are such 
that c((Ai,r)~-^,p~^) = (rii,r, ri2, p), that is, the definition (4.9) of the reduced variables 
holds. More precisely, (.np has the nonlocal expression 

^np = II f/(e(s, s), k{s, s'),p{s),p{s'),T{s),T{s'))dsds', (4.10) 

where U : S0{3) x (M'^)^ — ;■ M is a given function and 

k{s,s') ■.= -A^\s){r{s)-r{s')) G M^, 
as,s'):=AY\s)Aiis')G50i3). 

This general expression of inp allows the treatment of a potential depending on the Euclidean 
distance between two charges. Indeed, as we can see in Figure 1.1, the position of a charge is 
given by 

Ckikiis) = r{s) + Ai(s)T7fc^(s) + A2{s)Ai{s)r]i,^{s), 

and so the distance between two charges is found to be 

Ofci/c2"il"i2 \'^' ^ / — I ^k\k2 \Sj ^kik'z (^ ) I 

= |r(s) + Ai(s)r7fc^(s) + A2(s)Ai(s)r7,,^(s) -r(s') - Ai(s')r7„,(s') - A2(s')Ai(s')t7^,(s')| 
= \A^\s){r{s)-ris'))+r„^is) + AY\s)A2is)A,is)r„^is)- 

AY\s)A^is')rJ^^is') - A^\s)A2is')A^is')r,^^is')\ 
= \-Kis, s') + rjk, (s) + p(s)T7fc^ is) - as, s'Wrm («') - e(s, s)p{s')r,^^ {s') | . (4.11) 

This expression shows that the distance between two charges can be expressed solely in terms 
of the variables ^(s, s'), k(s, s'),p(s),p(s'). In addition, we note that Inp is invariant under 
the {Gi)l = 50(3)-action 

(Ai(s),r(s),A2(s)) ^ (/iAi(s),/ir(s),/iA2(s)/i-i), h £ 50(3). 

Hence, we are in the situation described in §3.4. 
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Reduced equations of motion. In order to write the reduced equations (3.5), we need 
to consider the affine representation 9 of the group G = T{I,SE{3)(§)SO{3)) on V* = 
J7i(/,se(3)(Dso(3)) X T{I,R^), given by 

Vi,r,A2)("i' r, ^2, p) = (Al, r, A2) (fii, r, ^2, p) + c(Ai, r, A2), 

where the first term denotes the representation of G on V* defined by 

(Al, r, A2)(ni, r, n2,p) = (Aifii, Air - Aifii X r, A2Ai(ni + ^2) - Aifii, Aip) . 

Thus, upon using the abstract system (3.5), the reduced Euler-Lagrange equations become 



5j 



6i ^ 5i 6£ _i _.6e 
+ CJ2 X - — h ri2 X — — + —p ' -p ^ = 

00^2 05 '2 op op 



{dt + ui x^— + {ds + ^1 y<^ 



(4.12) 



5-^ \ J 6T 6p 

9f +0^2 X h 9s + $12 X TT^+CJl X hfii X — — -p ^ =0. 

/ OU2 V / 0il2 OU2 0il2 op 

In these equations the functional partial derivatives are total derivatives. These equations 
could have been obtained also directly by using the explicit expressions of the variations for 
all variables, namely, 

(^wi = — Si + [wi,Si], ^7 = — * + c^i X * + 7X Si 

d 
SuJ2 = -WL^2 + Pi,a;2] + p2,wi] + p2,ii^2], 5p = pT.i -T.ip + p'E2 

50i = — Si + [f7i,Si], (^r = — * + f2i x* + rxSi 

OS OS 

5^2 = I-S2 + pi, (^2] + P2, J^l] + [S2, ^2], Sp = pX^i+ *. 

as 

where Si := A^^Mi, * := A^'^Sr, S2 := A^1A2 ^M2Ai. 

If we use the explicit expression of ^ = iioc + inp, with inp given in (4.10), this system 



Gay-Balmaz et al. 



Exact geometric theory of dendronized polymer dynamics 



24 



becomes 



dt + ujix ) \- (ds + ilix ) ^j^ + p X 



+ 7 X 



6uji 



<5ni 



5j 






^2 X -? h ^2 X ^^ 1 ^ P 

0U2 Olt2 op 

(^^(^'^'Ms)"' -P(s)"'^(s,s')j ds' 
("1^(5, s') X /^(s, s') + Z{s, s')] ds' = 0, 



Hoc ~r ^npj 

6V 
5p 



(4.13) 



9f + cji X j — ^ + f 9^ + fii X j 



d'y 



5T 



5p 



/( 



i{s,s'f^{s\s)-^^{s,s')]ds' = Q, 



dt+u:2X ] -j^ + (ds + fl2>i 



dn 

5koc 
5^0 



+ CJi X 



(5a;; 



+ fii X 



51 



loc 



6^2 



Pi^r'^is, s') + p{s)-'f^{s', s)) ds' = 0, 
dpi Op2 J 

where the functional derivatives of £ioc and inp are now usual ones and the term Z(s, s') G M'^ 
is the vector corresponding to 

^ f)TT f)TT 

Z{s, s') = ^{s, s')j^{s', s) - —{s, s')e{s, s') G 50(3). 

The expressions dU/d^, dU/dn, dU/dpi, and dU/dpi denote the partial derivative of f/ = 
U{^, K,pi,p2) as a function on 50(3) x M'^ x 5*0(3) x 50(3). We identify the cotangent space 
T2^SO{3) with r^50(3) by using the pairing 

(Va, C/a) := V • u, for ^ = A'^Va, u = -A-^Ua. (4.14) 

Using this identification, we have dU/d^ £ T^S0{3) and d/dpi S rp^50(3). 

To these equations one needs to add the kinematic equations (3.6) which in this case take 
the forms 

p = pu>2 — <Ji>lP + p^i 

{dt + a;ix)rii = dsUJi 

{dt + u;ix)T={ds + ftix)-f 

{dt + (^1 +a;2)x)i72 = {ds + rJix)cj2 
{dt + LL>ix)p = 'y. 

Note that the terms involving ^ and k do not enter the angular momentum for the second 
bouquet. This is because there is no term analogous to ^ involving A2. 

The Poisson brackets for the Hamiltonian evolution of two-bouquets may be obtained 
from the general formula given in §3.3. 



4.2 Computation of derivatives for the nonlocal term 

While equations (4.13) define the complete equations of motion for a two- level bouquet string, 
for practical applications it is useful to provide explicit formulas for the partial derivatives in 
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that equation and clarify their physical meaning. The terms 

Siloc Slloc 

are the angular and linear momenta in the coordinate frame associated with the filament. 
These momenta are connected to a;i and 7 by the geometry of the filament and are usually 
linear functions of Ui and -f. Similarly, 



are the elastic deformations associated with twisting and stretching of the filament. These 
quantities are usually related by the constitutive relations to the deformations themselves, 
which defines the physics of the base elastic rod. The terms involving the derivatives 

are the angular momentum and torque applied to the bouquet, respectively, and are related 
by physical and geometric considerations to (cjj, 7) for the first term and to (ilj, T) for the 
second term. In any case, the derivatives of the local terms should not cause any mathematical 
difficulty (although they may require quite detailed physical considerations). 

More care must be exercised in computing the derivatives with respect to the nonlocal 
variables {^,K,p). The derivatives of the local (elastic) Lagrangian 

op op 

denote the elastic deformations of second level of bouquets with respect to the first level. 
The p-dependence of the local Lagrangian may also come from geometric factors affecting 
moment of inertia of a two-bouquet, or total potential energy of the deformation. Without 
considering a particular molecule, these terms are difficult to identify. On the other hand, 
using the procedure outline below, these derivatives should be relatively easy to compute 
once the functional dependence of iioc on p is known. 

We therefore concentrate on the derivatives in the nonlocal term, which is both very 
instructive and can be applied to a general potential in view of concrete applications. 

We are now ready to compute the variations in the nonlocal terms. For convenience, we 
define the vector 

^kik2mim2\^i ^ ) 

= -k{s, s) + r7fc^ (s) + p{s)T]k^ (s) - <e(s, s')r]^^ {s') - C{s, s)p{s')r]^^ {s') (4.15) 

with 

(^kik2mim2\^ } ^ ) |^fciA;2»nim,2 V'^5 ^ / | • (_4.iDj 

In applications, the potential energy is a function of the distances dkik2mim2{s,s') between 
any two given charges. More precisely, if U{d) is such a potential (that could be electrostatic, 
screened electrostatic, Lennard-Jones, or any combinations of those or any other potentials), 
then the nonlocal part of the Lagrangian is given by 



"np — 7 ^ 

ki,k2,mi,m2 



ju{dk,k2m,m2{s,s'))\T{s)\\T{s')\dsds'. (4.17) 
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The reverse distance between the points s and s' is determined by the following formula: 

= (.''^{s, s')k{s, s') + r]^.^{s') + p{s')r},,^{s') - CHs, s')t]^^{s) - C^{s, s')p{s)r]^^{s) 

= -C^{s, s')dm^rn2kik2{s, s') (4.18) 

The derivatives with respect to S, and p are computed using the following lemma. 
Lemma 4.3 Let f : S0{3) — >■ M 6e the function defined by 

/(A) = l|b-Aa|2, 

where a, b G M'^ are constant vectors. Then the partial derivative of f , relative to the pairing 
(4.14), is given by 

-4 = (b X Aa^A = (d X ^a)^^, 
dA ^ ' ^ ' 

where d = b — via. 

Proof. Given Va G ^^50(3) and v = A~^Va, let C(e) be a smooth curve satisfying C(0) = A 
and C'{0) = Va. We have 



5/ ^ \ _ ^ 
dA'^^f-Jt 



f{C{e)) = - (b - Aa) • Vas. = - {A'^h - a) • (v x a) 

e=0 

(^-^b X a) • V. 



Therefore, from the definition of the duality pairing (4.14), we conclude that 

df 



A{A'^h X a)^ = (b X AayA, 



as desired. 



Using this formula with d = dk^i.^mim2{^i ^') ^^^ A = ^(s, s'), we deduce that the partial 
derivative of U relative to ^ is given by 

dU 

= E ^'^"^'^'^-^-f^f^^ (d,,,,^,^, X e (^7^,(^0 +p{s')Vm2is'))Y^, (4.19) 

, , "'kik2mim2\^: ^ ) ^ ' 

kik2mim2 i ^ i ^ \ ' 

where dk^k2mim2 ^-nd ^ are evaluated at {s,s'). Using (4.18), we get 
dU 

= Y^ fcifcamima ^ , M - T^d^imafcifca X T^ (^mi (^) + P(^)^m2 (g))) ^C^ 

= E ^y'^'^"^^";!f:f r^(-d^,^,.,., X {v^M+P{s)Vm2{s))Y. (4.20) 

kik2mim2 
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where dmjm2kik2 ^iid ^^^ are evaluated at (s,s'). Thus the vector Z{s,s') appearing in the 
equations of motion can be exphcitly computed using the formulas 



^^(■s,s)C (s,s) 



^ \(J'kik2mim2\^i ^ ) ) 






C(s>'s')-^(s',s) 

Terms involving derivatives with respect to p are computed similarly. We are not going to 
present them all. For example, 

^"^('')7w;Y^^'''')= Y1 fcifc2mim2 , (p-\s)dk,k2mim2is,s) XTj^^Gs))^ 

'^Py^J , , "'fcifc2mim2 v*) * 7 



kik2mim2 






5 Generalization to A^-bouquets 

We shall now turn our attention to the general case of A^-bouquets. For the generalization 
to A'^-bouquets, the position of the charge is given by 

Cfcifc2...fciv(s) = r{s) + Ai(s)77fc^(s) + A2(s)Ai(s)r/fc2(s) + ••• + AAr(s)...A2(s)Ai (5)77^^(5), 

where r{s) is the parametrization of the filament, tj^Xs) is a vector of constant length that 
determines the position of the the branches at level i, where ki are the labels of the branches 
at this level, all of this in the reference configuration. The proper rotation matrices Aj(s) give 
the orientation of the ith branch at parameter value s. Generalizing the case of 2-bouquets, 
the potential depends on the distances 

^kik2..kff,mim2...mN\^T^ ) • |'^fcifc2...fejv W/ Crrai-m2---'"ijv V'^ )\ 

between the charges (fei, k2, ■ ■ ■ , k]\f) and (mi, m2, • • • , wijv) on the filament at position s and 
s' , respectively. 

We have S'0(3)-invariance of the distances under the transformation 

(Ai(s), r(s), A2(s), . . . , Aiv(s)) ^ (Mi(s), hr{s), /iA2(s)/i"\ . . . , /iA^(s)/i"^) 

for any h G S0{3), which shows that SO{3) acts on the iterated semidirect product 

SE{3)®{SO{3){(i)---SO{3))---) with (A^ - 1) copies of 50(3). 

We apply the theory developed in §3.5 for d = F{I, SE{3)), G2 = . . . = G^ = T{I, 50(3)). 
Thus the Lagrangian variables are 

(Al, r, A2,. ..,A^)G^(/, 5^(3) (D (50(3) (D(...(D 50(3)))) 
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and the reduced variables are 



Wl 


= 


Ar'Ai 




fii = 


= Ar^A; 


7 


= 


A^^f 




r 


= A^V 


UJ2 


= 


A^U^UsAi 




^2 -- 


= A^iA^iA'2Ai 


W3 


= 


Ar^A^iA^^AsAaAi 




^3 -- 


= A^U^^AgU'gAzAl 


UJN 


^ 


A^^ • • • A]^UjvAjv_i • 


••Al 


^N = 


= A^^...AJj^A'j^An-i ■ 


P2 


= 


Ar^AsAi 




P 


= A^V. 


P3 


= 


A^U^UgAaAi 








PN 


= 


Aj;^..AAr...Ai 









■Al 



(5.1) 



Remark 5.1 In the general case discussed above, there will be indices associated with A^*" 
labeling the attachment point of each bouquet frame. Then, the reduced variables in (5.1) 
will acquire indices as well, for example 

p^i = A^^A^^Ai 

^^1^=2 =A^i(A^i)"^A§2A^iAi 

p5;'^-'^-^=Ar^..A5f-\..Ai. 

As we discussed earlier, this indexing poses no fundamental mathematical difficulties, but 
encumbers the notation in the resulting formulas. For the sake of simplicity and clarity of 
the exposition, we shall avoid this indexing in our further discussions. ♦ 

The material Lagrangian depends only on (Ai, Ai, r, r, A2, A2, . . . , A^v, A^r). As in the case 
of 2-bouquets the Lagrangian is the restriction of a parameter-dependent Lagrangian de- 
fined on TT{I, SE{3) (D (50(3) (S) (• • • (S) 50(3) •••))) x T{I, R^)^+^ for the zero value of the 
parameter in J^(/, M'^) "*"'^, that is, 

L{ao=o) = L{Ai,Ai,r, r, A2, A2, . . . , Aat, An). 

We are thus in the situation described by the theory developed in §3.4, provided L is (Gi)o- 
invariant. We have (Gi)q = 50(3) C T{I,SE{3)) as a direct verification using (4.8) shows, 
when extended to A'" copies of 50(3). 

The most general form of the symmetry-reduced Lagrangian for A^-bouquets is 

i =iioc{i^l,'r, ^2, • • • , ^N,P2, . . . PAT, r2i, r, ^2, . . • , ^N, p) 

+ 4p(a;i, 7, ^2, • • • , ^Af,P2, • --PN, rii, r, ri2, . . . , ^n, P, (Ai, r)), 

where iioc is explicitly given in terms of (0^1,7,0^2, • • • ,^n,P2, ■■ -Pn, ^1,^,^2, ■■ ■ ,^n,p) 
and inp has still a dependence on (Ai,r), where (Ai,r) are such that (5.1) holds. More 
precisely, inp has the nonlocal expression 



^np 



J J U{as, s'),Kis, S'),P2{S), . . . ,PN{s),P2is'), . . . ,PN{s'),T{s),T{s'))ds ds'. 



where U : 50(3) x (R3)^+i ^ R is a given function and k{s, s') £ R^ ^(s, s') £ 50(3) have 
the same formulas as for the 2-bouquets. 
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Given a Lagrangian i, formula (3.18) yields the following motion equations for A^-bouquets 
comprising {N + 1) equations: 

, Si ^ 6i ,^ ^ , 6i ^^ 6i 
(dt + uix)- ^l^^k X T \-{ds + flix) —— + 2^rifc X 



5u>i 
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5uk 
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61 ^61 51 ^^ 

+px— +rx^:^+7x— + 2^ 



5p 
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6i 



k=2 



Sftk 



.1 se 



k^z 



{dt + uix) — + a, + fiix — - — = 
07 oi op 



, 6i '^ 6i ^ 51 

[dt+U,x) — + ^UkX—+ }_^ U^kX 



k=l 



k=i+l 
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5uk 
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^r, r^ , Se ^^ 5i s^ ^ 5i 



5ni 
* fc=i 
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^^^"-^^'^,'=° 
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+ (5s + ri/vT X 



^cu 
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5^ 



+ > , fifc X 



Af 



A:=l 



(5n 



N 



■Pn 



5p 



N 



(5.2) 



where all the partial functional derivatives are total derivatives. Note that x^(s) S T*,.S0{3) 

and hence Pi{s)~^-^{s) G so(3)* = M'^. As in the case of 2-bouquets these equations can be 
made explicit and one gets an analogue of the system (4.13). 



6 The equations of motion in conservative form 

In order to elucidate the physical meaning of dynamic equations (3.7), we shall recast them in 
conservative form. In principle, the form of these equations express the laws of conservation 
of linear and angular momenta at the every branch level. However, the direct formulation 
in terms of these physical quantities leads to very cumbersome expressions, as the linear 
and angular motion of each branch contributes to the motion of all branches on the next 
levels, and the conservation laws for different bouquet levels are thus highly intertwined. It 
is preferable to use the geometric formulation and express the conservation laws as coadjoint 
motion on the corresponding semidirect product groups. In that case, the formulas become 
highly compact and yet their physical meaning remains clear. This geometric approach also 
lets us formulate the sequential conservation laws in §6.3 below. 
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6.1 Two-bouquets 

Using (3.7) and (4.7) , we see that the equations of motion (4.13) take the form 

a^_ \ /^ M _5r\ ( d__ \ ( M^ 51 M^ 

dt ^Vi,7,-2)J [su,'6-r' 6U2) ^ [ds ^Vi.r,n.)J [sfti' 6T' 6^2 

— xp,— ,0 +U V-:rP '0'^ V 

op op J \ op op op 

Therefore, we can write these equations in conservative form by using (3.10). We get 

d f 6i 6i 6i\ d ( bl 6i 5i 

at Vi,-,A2)- [su,'d-f' 6u;2) ^ 9>.^Vi,r-,A2)- [sQi'dT' 6^2 



Ad{Ai,r-,A2)-i 



6e 



X P 



5i 



6p 5p 



+ U 



-1 



5p 5p 



1 r. 1^^\ 



(6.1) 



(6.2) 



6.2 A^-bouquets 

For A^ bouquets, we obtain from (3.19) the following conservative form 



d 



AdT. 



6£ 6£ 



se 



^^ — (Ai,r,A2,...,A^)-i \^5^^' 5^'-' §^^J + 5^ Ad^Ai,rA2,...,A^)- 
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6ni 6T' ' dn 
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Ad 



(Ai,r,A2,...,Ajv)~i 
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,1 d£ 
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TV 



_id£ -^ f 6£ _i _i S£ 



_,5£ ^ 



5pi 



k=i+l 



Spk^'' ^^ 6pk 



\ -Pn 



6£ 



dp 



'N 



(6.3) 



6.3 Sequential conservative forms 

It is interesting to cast the equations in conservative form not just for the strand as a whole, 
but for each particular bouquet and all the corresponding sequential bouquets. As an everyday 
analogy, consider a tree with many branches swinging because some force shakes the trunk. 
The conservative form for the whole tree sums all the momenta for each branch and proves 
that all those momenta would sum to zero, were it not for the external force acting on the 
trunk. 

Let us now introduce an imaginary cut to a tree branch, and try to apply the corresponding 
conservation laws. All the forces (and torques) starting from the branch and propagating 
above to smaller branches must sum to zero in a conservation law. The external force will 
be now applied by the tree to the base of this branch. We can now continue to higher and 
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higher branches and, as long as the tree structure is completely sequential, we should find an 
exact conservation law for all the branches at the level n = 1,2,... and higher. 

Comparing the conservation laws for A^-bouquets (6.3) with the equations of motion (5.2) 
for i^-bouquets, we find a new conservation law for the first K branches in the A^-bouquets 

6£ 6£ 

— X p, 

(5p Sp 

N 



Ad^Ai,r,A2,...,AK)-M yrXP'TT'O'-'O 



y j=K+i ■' ■> 

MAi.r.A, A^l-i (^1>0>^2, • • • ,5'i^) , 



where 



N 



with the convention that pi is absent. 

More generally, the conservation laws for the branches J, . . . ,K, l<J<K<N,m the 
A^-bouquet is 

61 



Ad^A„...,AK)-i •••' E ^jX77r + "^x 



j=K+l ■> J 

( '\^ ^^ r^ ^^ 

(Aj,...,Ak) M ' Z^ 3 SuJi ^ S^i 

\ i=i 
For example, taking J = K = N = 2 for 2-bouquets yields 

V ' 

1st bouquet acting on 2nd bouquet 

7 Conclusion 

As with many other complex biological and organic molecules, the final conformation of the 
compound polymers considered here cannot be determined only on the basis of the forces 
acting on the polymer. The choice of the dynamical path in configuration space taken by the 
molecule during conformation is also of crucial importance. In this paper, we have outlined 
a geometrically exact theory of dendritic polymer dynamics based on symmetry reduction 
of Hamilton's principle defined on the tangent space of the iterated semidirect-product Lie 
group 

5^(3) (S) (50(3) (S) (S0(3) ©{■■■® 50(3)) • • • )) 

for Lagrangians that are invariant under changes of orientation SO{2>) in the first Lie group, 
SE{3) ~ 50(3) (DM^. As in the example of shaking the trunk of a tree with many branches. 
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the overall orientation and position of the tree may be immaterial, but the relative orientations 
of its successive branches are important to the time-dependent oscillations and deformations 
of the tree away from its reference configuration. The motivation for the iterated semidirect 
product action was illustrated by considering the case of coupled rigid bodies, and was de- 
veloped into a general theory of dendritic polymers with many branches and types of rigid 
charge bouquets. However, specific configurations and individual molecular parameters, such 
as moments of inertia and length-scales, were not considered. Although the approach was 
general, the calculations were performed explicitly enough for real applications of evolutionary 
modeling of dendritic polymers to be made, once their reference configurations and molecular 
parameters were specified. We concentrated mainly on two-level bouquets, because this case 
illustrates most of the mathematical concepts and is the necessary first step in applications. 
However, formulas for A^-level dendrimers were also derived and were discussed in enough de- 
tail to make their general pattern evident for understanding the scope of further applications. 
Both local and nonlocal interactions were included in the theory, in order to account for the 
screened electrostatic interactions that occur among different branches of the dendronized 
polymers. Finally, the theory was expressed in conservative form, so that forces and torques 
among the branches could be explicitly identified. 

While our theory is very general and is applicable to a wide range of elastic and nonlocal 
forces, the major drawback at this stage is that it does not consider dissipation due to possible 
presence of ambient fiuid. Indeed, the most natural media for polymers is a liquid solution, 
where the forces of friction dominate the dynamics. Our theory does not consider the forces 
of friction and dissipation and thus is not readily applicable to polymers in liquids. It is not 
clear at this stage how to derive geometrically exact analytical expressions for friction forces 
acting on the dendrimers in an arbitrary conformation. Presumably, further progress in that 
field will need to include some insightful, yet ad-hoc, physical approximations for frictional 
forces as a function of local velocities in the semidirect product groups considered here. We 
did not dare to venture into this territory here and only considered the case when the friction 
effect due to the ambient media can be neglected. 



8 Acknowledgements 



We thank M. Bruveris, D.C.P. Ellis, and J.E. Marsden for helpful discussions of this work as it 
was being developed. The work of FGB was partially supported by a Swiss NSF postdoctoral 
fellowship. DDH and VP were partially supported by NSF grants NSF-DMS-05377891 and 
NSF-DMS-09087551. The work of DDH was also partially supported by the Royal Society of 
London Wolfson Research Merit Award. TR was partially supported by a Swiss NSF grant 
and the program Symplectic and Contact Geometry and Topology at MSRI Berkeley. 

References 

[1] D. A. Tomalia and P. M. Kirchhoff. Rod-shaped dendrimers, 694. Midland, MI, USA: 
The Dow Chemical Company, US patent 4,694,064, 1987. 

[2] D. A. Tomalia, A. M. Naylor, and W. A. Goddard III. Starburst dendrimers: Molecular- 
level control of size, shape, surface chemistry, topology and fiexibility from atoms to 
macroscopic matter. Angew. Chem., Int. Ed., 29:138-175, 1990. 

[3] A. D. Schliiter. Dendrimers with Polymeric Core: Towards Nanocy Under s, pages 165- 
191. Topics in Current Chemistry. Springer (Berlin), 1998. 



Gay-Balniaz et al. Exact geometric theory of dendronized polymer dynamics 33 

[4] H. Frauenrauth. Dendronized polymers - building a new bridge from molecules to 
nanoscopic objects. Progress in Polymer Science, 30:325-384, 2005. 

[5] B. M. Rosen, C. J. Wilson, D. A. Wilson, M. Peterca, M. R. Imam, and V. Percec. 
Dendron-mediated self-assembly, disassembly, and self-organization of complex systems. 
Chem. Rev., 109:6275-6540, 2009. 

[6] D. J. Dichmann, Y. Li, and J. H. Haddocks. Hamiltonian Formulation and Symmetries 
in Rod Mechanics, volume 82: Mathematical Approaches to Biomolecular Structure and 
Dynamics, chapter II, pages 71-114. Springer IMA, New York, 1992. 

[7] D. K. Christopoulos, A. F. Terzis, A. G. Vanakaras, and D. J. Photinos. Helix formation 
in linear achiral dendronized polymers: A computer simulation study. J. Chem. Phys., 
125:204907, 2006. 

[8] Y. Ding, H. C. Ottinger, A. D. Schluter, and M. Kroger. From atomistic simulation to 
the dynamics, structure and helical network formation of dendronized polymers: The 
janus chain model. J. Chem. Phys., 127:094904, 2007. 

[9] R. Connolly, G. Bellesia, E. G. Timoshenko, Y. A. Kuznetsov, S. Elli, and F. Ganazzoli. 
intrinsic and topological stiffness in branched polymers. Macromolecules, 38:5288-5299, 
2005. 

[10] D. D. Holm, J. E. Marsden, and T. Ratiu. The Euler-Poincare equations and semidirect 
products with applications to continuum theories. Adv. Math., 137:1-81, 1998. 

[11] F. Gay-Balmaz and T. S. Ratiu. The geometric structure of complex fluids. Adv. in 
Appl. Math., 42:176-275, 2009. 

[12] J. C. Simo, J. E. Marsden, and P. S. Krishnaprasad. The Hamiltonian structure of 
nonlinear elasticity: The material and convective representations of solids, rods, and 
plates. Arch. Rat. Mech. Anal, 104:125-183, 1988. 

[13] D. D. Holm and V. Putkaradze. Nonlocal orientation-dependent dynamics of charged 
strands and ribbons. C R. Acad. Sci. Paris, Sir. I: Mathematique, 347:1093-1098, 2009. 

[14] D. Ellis, D. D. Holm, F. Gay-Balmaz, V. Putkaradze, and T. Ratiu. Geometric mechanics 
of flexible strands of charged molecules. Arch. Rat. Mech. Anal., To appear, 2009. 

[15] F. Gay-Balmaz, D. D. Holm, and T. S. Ratiu. Geometric dynamics of optimization. 
preprint, 2010. 

[16] J. E. Marsden and T. S. Ratiu. Introduction to Mechanics and Symmetry: A Basic 
Exposition of Classical Mechanical Systems (Texts in Applied Mathematics). Springer, 
2nd edition, 2002. 

[17] R. Grossman, P. S. Krishnaprasad, and J. E. Marsden. The dynamics of two coupled 
rigid bodies. Phys. Lett. A, 133:134-139, 1988. 



